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1 .  Introduction 

In  a  recent  paper  by  G.  Dantzig  [1],  an  algorithm  was  introduced  to 
obtain  solutions  for  certain  quadratic  programming  problems  by  a  modified 
simplex  program.  This  programming  algorithm  follows  closely  the  methods 
employed  by  Wolfe  and  Markowitz  [2].  In  this  introduction,  we  shall  out¬ 
line  Dantzig' s  theorems  for  those  unfamiliar  with  this  method,  and  also 
develop  the  notation  used  in  this  paper. 


1.1.  The  Criterion  Function 

We  shall  formulate  the  problem  as  the  constrained  maximization  of  a 
convex  quadratic  function  in  x.  The  criterion  function  (in  vector/matrix 
notation)  will  be 

(1.1)  n(x)  -  xS)  -  x*Ax, 

where 

x  >  0, 

t 

x  Ax  >  0  for  all  x,  and  A  is  a  symmetric  matrix. 


This  research  was  supported  by  Varian  Associates,  Inc.,  Palo  Alto, 
California.  Publication  was  supported  by  Office  of  Naval  Research  Contract 
Nonr-225(50)  at  Stanford  University.  Reproduction  in  whole  or  in  part  is 
permitted  for  any  purpose  of  the  United  States  Government. 

I  wish  to  thank  Mr.  G.  Kuykendall  and  Mrs.  L.  Owens  of  I.B.M. ,  San  Jose, 
California,  for  their  help  in  introducing  the  author  to  1401  Fortran  and 
in  testing  this  program. 


ggggggaga  bwrww 


Some  of  the  columns  of  the  A  matrix  may  have  zero  elements.  The  activ¬ 
ities  associated  with  these  zero  rows  and  columns  will  be  linear  activities, 
not  quadratic  activities  Thus,  by  suitable  construction  of  the  A  matrix, 
mixed  quadratic/linear  solutions  may  be  obtained. 

1  2 .  The  Constraint  Function 

Tn  addition  to  the  positivity  requirement  for  x,  we  require  that 
the  following  equation  (in  vector/raatrix  notation)  be  satisfied: 

(1.2)  Cx  =  d  . 

1  3  A  Note  on  Slack  Variables 

The  vector  x  is  defined  so  that  it  contains  the  required  slack 
components  to  make  (1.2)  hold  as  an  equality. 

1  u  Tne  Kuhn-Tucker  (3]  Optimality  Conditions 

A  solution  to  (1  1)  and  (1.2),  toy  the  vector  x*.  can  be  segregated 
into  two  subsets.  Let  X°  be  the  subset  of  components  of  x  which  are 
zero,  and  X+  the  subset  of  components  of  x  which  are  positive.  Tne 
solution.  x»,  is  the  optimal  solution  (leads  to  a  maximum  jt(x)  under 

constraints,  (i  2))  if 


(1  3) 

b  -  Ax* 

-  C“h*  =  0  for  x*  e  X* 

(1  “) 

b  -  Ax* 

•  c\*  <  0  for  x*  t  Xu  , 

• 

and 

(1-5) 

Cx*  *  d 

The  vector  k*  is  defined  as  the  Lagrangian  multiplier  vector. 

1  5 •  The  Complementary  Vector 

Suppcse  we  define  a  new  vector  u*  such  that 
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(1.6) 

and  such  that 


u#tx  =0  ,  u*  >  0  , 


(1<7)  b  +  u*  -  Ax*  -  c\*  =  0  for  x*  e  X'  . 

We  can  now  write  ( 1 • 3 )  and  (l-M  as 

8)  b  +  u*  -  Ax*  -  CtA*  =0  for  all  x  >  0  . 

It  can  be  seen  then  that  the  vector  u  is  a  special  kind  of  slack 
variable  which  permits  the  Kuhn-Tucker  inequality  (l.U),  together  with  U-3)> 

to  be  written  as  an  equality  (1*6). 

1  6  Similarity  Between  These  Equations  and  Linear  Programming  Equation^ 

The  similarity  between  the  satisfaction  of  equations  (1.5)  w»d  (1.8) 
and  simplex  method  linear  programming  can  be  exploited  if  satisfaction  of 
(1.6)  can  be  ensured  and  some  method  of  determination  can  be  made  of  the 
efficient  order  in  which  to  bring  in  new  components  of  x  (to  take  the 
place  of  the  linear  criterion  function  coefficients  of  linear  programming). 
Dantzig  has  supplied  such  a  selection  criterion.  A  feasible,  non-optimal 
solution,  say  x°  (which  satisfies  (1-5)  but  does  not  satisfy  (1-3)  or 
(l.t*)),  will  be  characterized  by  a  u  vector  which  contains  negative 
components.  The  component  of  x°  corresponding  to  the  most  negative 
component  of  u°  is  the  activity  to  be  introduced  next  into  the  basis 
(X+).  Dantzig  has  shown  that  such  a  selection  rule  will  lead  to  a 
mono tonic ally  increasing  n(x). 

2.  Setting  Up  the  Solution  Matrix 

2,i,  The  Complete  Form  of  the  Solution  Matrix 

Figure  2.1  shows  the  solution  matrix  for  the  quadratic  criterion 
function  (1.1),  subject  to  the  explicit  constraints.  Note  the  requirement 
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so  that  Fig. 2.1.1  may  be  represented  as 


One  practical  disadvantage  of  the  simplex  quadratic  programming 
algorithm  is  the  very  large  size  of  the  solution  matrix.  With  m  activities 
and  n  constraints,  the  solution  matrix  is  a  (m+n)  x  (2m+n)  matrix,  where 
m  >  n.  This  large  matrix  limits  the  application  of  the  method  to  either  the 
solution  of  small  problems  or  the  use  of  large  computers. 

We  note  that  n  of  the  &  activities  are  slack  activities,  the  rest 
being  productive  activities.  We  may  partition  the  x  vector  ao  follows: 

(2.2)  X  =  <  xp;*s  > 

where  x  contains  the  components  of  x  which  represent  slack  activities, 

x  contains  the  components  (m-n  in  number)  of  x  which  represent 
P 

productive  activities. 

In  the  same  way,  we  write 


Furthermore,  we  will  assume  that  the  criterion  function  is  restricted  only 
tc  the  productive  activities;  that  is,  the  slack  activities  are  not  costly 
or  productive.  Thus,  we  have 


(2.M 


b  =  <  b  >  , 

r 


where  b  contains  only  the  coefficients 

P 

activities. 

Also,  we  have 


(2.5) 


A  = 


'_Ap  :_0_; 

.  0  ;  0  . 


corresponding  to  the  productive 


where  A  is  the  matrix  whose  elements  are  associated  only  with  the  produc- 
P 

tive  activities,  and  r  p 


(2.6) 
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where  C  is  the  matrix  associated  with  the  productive  activities,  and 
P 

I  is  the  unit  matrix  associated  with  the  slack  activities. 

The  whole  solution  matrix  of  (2.1)  may  he  rewritten  as 


(2.7) 


b  1  0  1  d 

P  1 

AP  j  *  !  CP 

X 

p 

0  10*1 

X 

s 

■ 

_ 

ts. 

T9. 

u 

p 

*  \  -  i  * 

u 

s 

<=;  | ! :  f 

X 

We  note  from  the  last  two  rows  of  (2-7)  that  ug  and  X  are  related  in  a 
very  useful  way.  An  increase  in  a  component  of  u^  implies  an  equal  increase 
in  a  corresponding  component  of  X.  and  vice  versa.  Ihius,  for  any  component 
in  u  and  X,  we  may  write 

S 


(2-8) 


uj +  6j  =  +  V  J  =  1,2'*“'n  ’ 


where  &  is  some  change  to  u,  or  X,.  For  all  practical  purposes,  any 
j  J  J 

selection  criterion  bused  (in  part)  on  u  could  be  based  also  on  X.  The 

5 

solution  matrix  (2.7)  may  then  be  rewritten  in  a  reduced  form  We  have 


(2  9) 


b  1  d 

P  1 

A  i  C 

P  ,  P 

X 

P 

0  1  I 

— f - 

! 

W  1 
X  I 
_ L 

-I  *  0 

U 

_J2 _ 

C*  ’  0 

p  >  * 

X 

Where  we  have  m  activities  and  n  constraints,  the  general  solution  matrix 
of  dimension  (m+n)  x  (2m+n)  may  be  reduced  to  a  m  x  2m  matrix.  Such  a 
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reduction  leads  to  a  marked  saving  in  computer  memory  elements. 


3.  Ihe  Computer  Program 

3 . 1  General  Information 

The  computer  program  is  written  in  Fortran  language  fox  the  l.B.M. 

1401  computer  f4V  Ihe  1401  computer  is  a  medium-sited  laaehioe  which  is 
quite  cottraoaly  fonid  if!  use  as  an  accounting  maetioa  in  many  thcusaruis  of 
business  firms  ia  the  Uhitad  States  and  abroad.  £iaee  quadratic  program¬ 
ming  is  a  vaiaat>*s  aid  tc  profit  anaujsis  it  fcuaiuess,  we  feel  that  thj.s 
program  will  find  oar y  uses  is  these  firms.  Approninateijr  8,030  memory 
positions  ere  required  tu  store  the  progras  and  the  necessary  fixed  routine 
instructions  for  the  computer.  Thus,  if  12,000  memory  petitions  are  avail¬ 
able,  sufficient  memory  positions  remain  to  store  a  15  1  3®  solution  matrix 
with  an  accuracy  of  sis  si gx<if 5 cart  figures  (this  includes  space  for  the 
vectors  x  .  x  ,  u  ,  t  ,  and  d>.  This  solution  capacity  should  be 

p  5  p  p 

adequate  for  most  business  problems.  One  solution  cycle  requires  approx¬ 
imately  *0  seconds  of  machine  time  and .  o?;  the  average,  acout  10  cycles 
are  required  to  find  an  optimal  solution.  Thus,  a  solution  may  be  obtained, 
or  the  average,  in  5  minutes  of  computer  time.  Pc.  magnetic  tape  units  are 
required,  the  problem  being  solved  completely  in  the  magnetic  core  storage 
areas. 

3.2  Description  of  the  Block  Diagram 

It  is  necessary  to  make  the  computer  Fortran  program  notation 
different  from  the  descriptive  notation  used  in  the  previous  sections. 

The  x  vector,  u  vector,  and  X.  vector  are  known  to  the  computer  as  a 
single  X  vector.  The  assigned  index  numbers  are  the  only  clues  to 
identification  of  the  x,  u,  and  \  vectors  within  the  X  vector  Also,  the 
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c  and  d  vectors  are  known  to  the  computer  as  a  single  C  vector,  from 
which  the  c  and  d  vectors  may  be  identified  only  by  the  assigned  index 
numbers.  The  solution  matrix  is  known  to  the  computer  as  the  A  matrix. 

Hie  submatrices  described  in  previous  sections  are  only  identified  by  the 
assigned  index  numbers 

With  these  definitions  in  mind,  Fig.  3-2.1  may  be  examined  Figure 
3-2.1  is  a  simplified  block  diagram  of  the  modified  Dantzig  simplex  method. 
Each  block  will  be  described  in  detail  below 

Block  1 

During  this  operation  the  computer  accepts  input  information  from 
cards  and/or  tape.  Tills  information  is  placed  in  the  reserved  memory 
positions  and  the  assembly  of  the  A  matrix  is  completed  In  addition,  the 
C  vector  is  also  assembled  in  the  appropriate  memory  locations 

Block  2 

In  Block  2,  the  A  matrix  and  the  C  vector  ure  used  to  determine  the 
initial  value  of  the  X  vector.  This  initial  solution  is  farmed  by  running 
the  slack  activities  at  a  level  which  consumes  all  the  available  -escurces 
Furthermore,  the  complementary  productive  activities  are  a- so  introduced 
(since  the  productive  activities  ure  all  ct  zero)  so  that  the  C  vector  is 
completely  satisfied.  Hiis  initial  solution,  although  a  feasible  solution, 
is  far  from  optimal  since  all  the  resources  are  "thrown  away"  by  the  t».ack 
activities. 

Block  3 

The  printer  now  prints  out  the  initial  solution.  First,  the  C  vector 
is  printed  out.  Then  the  A  matrix  and  X  vector  (in  the  far  right  column) 
are  printed  out  row  after  row,  until  the  entire  matrix  and  vector  are 
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complete.  This  matrix  may  then  be  examined  to  assure  that  the  data  fed 
into  the  computer  has  been  stored  in  the  proper  positions. 

Block  i 

If  the  value  of  this  least-valued  complementary  activity  is  positive 
or  zero,  the  optimality  theorems  indicate  that  the  present  solution  is  the 
optimal  solution  and  the  computer  passes  on  to  Block  17.  Bovever,  this  is 
rarely  the  case,  especially  on  the  first  solution.  Ordinarily,  the  least¬ 
valued  complementary  activity  is  negative.  If  it  is  negative,  the  computer 
passes  on  to  Block  6. 

Block  6 

In  Block  6.  the  (most  profitable)  activity  is  increased  from  zero 
until  the  first  of  the  S  vector  constraints  coefficient  becomes  binding. 
Which  constraint  is  the  first  to  kind, is  determined  by  dividing  the  K  row 
of  A  matrix  elements,  one  at  a  ttse,  into  the  appropriate  C  recto*  coef¬ 
ficients.  Note  that  only  certain  C  creator  coefficients  (columns)  are  used 
during  this  process.  These  certain  coefficients  are  ia  the  coluer.s  asso¬ 
ciated  with  the  activities  which  are  restricted  to  sero  or  posit! ve  values 
orily;  that  is,  only  the  productive  activities  end  tfce  slack  activities. 

Tii us,  only  the  (?  vector  coefficient  associated  eith  a  sign-restricted 
activity  can  possibly  bind  the  incoming  activity.  the  complementary 
activities  and  the  Lagrangiah  ttultipliers  are  not  sign-restricted  for 
feasible  solutions.  One  exception  exists;  the  C  vector  coefficient  asso¬ 
ciated  with  the  complementary  activity  for  the  ktli  activity  is  sign- 
restricted  during  this  solution  cycle.  It  nay  not  be  permitted  to  go 
negative  at  this  time.  This  is  because  it  is  very  desirable  for  the 
incoming  activity,  I(k),  to  be  bound  by  the  C  vector  coefficient  which 
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formerly  bound  its  own  complementary  activity.  The  reasons  behind  this 
statement  vill  become  more  apparent  later  in  Block  13. 

Block  7 

Since  the  binding  C  vector  constraint,  L,  formerly  was  the  bound  for 
some  other  activity  that  was  formerly  in  the  basis,  this  other  activity 
must  be  identified  and  made  to  drop  out  of  the  new  basis  (go  to  zero) 
following  the  introduction  of  the  (most  profitable)  activity.  Thus 
activity  Is  found  by  finding  the  index  number  of  the  row  with  a  1 .  in 
column  L.  which  also  has  a  positive,  non-zero  X  vector  component.  This 
index  number  i  s  recorded  as  KA. 

Block  8 

Now  that  K,  L,  and  KA  are  known,  the  A  matrix  and  the  C  vector  are 
modified  according  to  the  simplex  rules.  Following  this  operation,  the 
Ktf  row  will  have  0.  everywhere  except  in  the  Lt‘  column  where  the  number 
1.  will  appear.  The  KAt!‘  row,  which  formerly  had  0. 's  everywhere  but  the 
L“”  column  (where  it  had  a  i.),  will  now  generally  contain  coefficients  of 
various  values.  Everywhere  else  the  A  matrix  coefficients  will  be  modified 
by  the  simplex  process  to  recognize  the  existence  of  a  new  basis  (contain¬ 
ing  the  K"  ,  instead  of  the  KA"  \  activity). 

Block  j 

At  this  stage,  the  new  X  vector  is  assembled.  The  incoming  activity, 
X(K),  will  become  positive  and  equal  to  the  new  C(K),  while  the  outgoing 
activity,  X(KA) ,  will  become  0.  .  Since  the  C  vector  has  teen  modified, 
the  other  activity  levels  will  change  somewhat. 

Block  10 

The  printer  will  print  the  sequence  number  of  the  new  solution. 
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Block  11 


The  printer  will  print  the  index  numbers  of  the  incoming  and  outgoing 
activities  and  the  binding  constraint. 

Block  12 

The  printer  will  cow  print  the  new  C  vector.  Following  this,  the  A 
matrix  and  1  vector  components  will  be  printed  out,  row  by  row,  until  the 
entire  solution  has  been  printed. 

Block  13 

The  conyplementari ty  check  is  now  made.  The  complementary  activity 
of  tee  incoming  activity  is  checked  to  see  if  it  is  0.  .  If  it  is,  this 
can  happen  in  the  following  way:  Suppose  that  the  binding  C  vector  coef¬ 
ficient  was  the  former  binding  coefficient  for  the  complementary  variable 
of  the  incoming  activity*  Xr.  this  case,  as  the  incoming  activity  comes  in, 
its  own  complementary  activity  goes  to  0.  .  Most  of  the  time  the  comple¬ 
mentary  activity  falls  fcut  does  not  go  to  xero  (this  is  because  some 
activity  other  than  the  -canyieoentary  one  is  dropped).  If  the  comple¬ 
mentary  activity  associated  with  the  in’coming  activity  becomes  zero,  the 
complementarity  test  passes.  If  ret,  the  test  fails.  The  optimality 
theorems  indicate  that  if  the  test  passes,  the  solution  Is  feasible.  If 
the  test  fails,  the  solution  is  ncr.-feasitle . 

if  the  complementarity  test  passes,  the  computer  moves  on  to  Block 
Ik.  If  the  test  fails,  the  computer  *oves  on  to  Block  15. 

Block  Ik 

She  computer  enters  this  block  only  when  the  complementarity  test 
passes.  The  printer  will  print  out  a  notice  that  the  solution  (printed, 
formerly  in  Block  12)  is  a  feasible  solution. 
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The  computer  now  returns  to  Block  4  and  determines  the  conditions 
for  another  solution. 

Block  15 

The  computer  enters  this  operation  only  when  the  '•omplementarity 
test  fails  and  the  solution  is  non-feasible.  According  to  the  optimality 
theorems,  the  complementary  activity  associated  with  the  binding  C  vector 
constraint  coefficient  must  be  brought  into  the  basis.  Therefore,  K  is  set 
equal  to  the  index  of  this  activity  and  the  computer  attempts  to  p-ass  the 
c cr.pl ementarl ty  test  again  by  going  back  to  Block  6,  with  this  new  »'alue 
of  K. 

Block  16 

This  block  is  entered  only  if  the  optimality  test  of  Block  5  is  passed. 
In  this  block  the  printer  prints  out  a  message  indicating  that  the  solution 
(printed  last)  is  the  optimum  solution.  Following  this  message  the  computer 
stops  and  the  problem  is  solved. 

3.  3.  Hoe  Fortran  Statement  of  the  Solution  Algorithm 

Appendix  I  is  a  reproduction  of  the  Fortran  statement  list  for  the 
solution  algorithm.  Tne  essential  heart  of  the  algorithm  starts  at  state¬ 
ment  sequence  number  30.  The  operations  preceding  this  point  are  the 
instructions  required  to  read  in  the  initial  data  and  form  the  initial 
solution  matrix.  These  instructions  c.an  sometimes  be  modified  to  suit  the 
particular  application  of  the  algorithm  as  a  component  of  a  .larger  program. 
In  this  case,  these  statements  have  been  made  general  enough  to  permit  the 
program  to  be  used  for  most  applications. 
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3.4.  The  Layout  of  Input  Data  Cards. 

In  the  assignment  of  coefficients  to  be  punched  in  the  input  data 
cards,  some  arbitrary,  but  prescribed,  index  order  of  the  constraints  and 
activities  is  assumed.  This  index  order,  once  set,  is  adhered  to  through¬ 
out  the  assignment. 

3.4.I.  The  Index  Carl 

The  index  information  is  placed  on  one  card  with  the  following 

format : 

Card  No.  Column  No.'s  Coefficient 

1  1-2  NT  (the  number  of 

productive 

activities) 

1  3-4  MT  (the  number  of 

constraints) 

The  following  equations  must  be  satisfied  when  selecting  NT  and  MT: 

i  <  NT  <  10 

1  <  MT  <  10 

2  <  NT  +  MT  <  15 


3.4.2.  Hie  C  Vector  Cards 

Hie  C  vector  is  placed  on  two  cards  with  the  following  format: 


Card  No. 

Column  No.'s 

Coefficient 

2 

1-8 

C(l) 

2 

9-16 

C(2) 

2 

17-24 

C(3) 

2 

25-32 

C(4) 

2 

33-40 

C(5) 

2 

41-48 

C(6) 

2 

49-56 

C(7) 

2 

57-64 

C(8) 

13 

3-«.3. 
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Card.  No. 

Column  No.'s 

Coefficient 

2 

65-72 

C(9) 

2 

73-80 

C(10) 

3 

1-0 

C(ll) 

3 

.  9-16 

C(12) 

^3 

17. Pa 

C(l3) 

3 

25-3? 

c(lt) 

3 

3  3- to 

c(l5) 

ed  components  of  C  are  left  blank. 

The  Quadratic  Coefficient  Cards 

We  ha,?e  UT  cards  of 

the  following 

format! 

Card  No. 

Column  No.  s 

Coefficient 

L 

1-8 

A(  1 » 1 ) 

u 

9 -*6 

A(1.2) 

1* 

17-2U 

A<  1 . 3 ) 

a 

25-32 

A(l.t) 

k 

33-to 

a  v.5) 

a 

1*1  -  i*3 

a(  I » 6 ) 

a 

ay- 56 

A(i,7) 

h 

57- 6a 

A<1,3) 

a 

65-72 

A (1 ,9) 

U 

73-80 

A(  .  .10) 

.e  above  sequence  fs 

repeated  for  NX  cards,  the  la: 

3  ♦  HI 

-73.80 

A( NT. 10) 

unused  A  eiemerts  are  left  blank. 

The  Constraint  Coefficient  Cards 

We  have  MT  cards  of 

the  following 

format. 

Card  No. 

Co 1 umn  No.'s 

Coeffi clent 

3  ♦  NT  ♦  1 

1-8 

c(i,i) 

3  ♦  NT  ♦  1 

9-i6 

C(2,l) 

3  ♦  NT  ♦  1 

17-2U 

C(3,l) 

la 


f  which 


Card  No. 

Column  No.'s 

Coefficient 

3  +  NT  +  1 

25-32 

C(4,l) 

3  +  NT  +  1 

33-40 

0(5,1) 

3  +  NT  +  1 

41-48 

c(6,i) 

3  +  NT  +  1 

49-56 

c(7,i) 

3  +  NT  +  1 

57-6* 

0(8,1) 

3  NX  ♦  1 

65-72 

0(9,1) 

3  ♦  NT  ♦  1 

73-80 

0(10. 1) 

repeated  for  MI 

cards,  the  last  of  which  is... 

3  ♦  NT  ♦  MT 

i  £ 

C(l.MT) 

3  ♦  KT  ♦  MP 

9-16 

C(2,MI) 

. . .  end  so  ot  until  . . 

• 

3  ♦  NT  ♦  MP 

73-30 

0(10, MT) 

All  the  unused  elements  of 

C  are  left  blank 

>• 

3.4.5.  Special  Instructions  on  Punching  These  Cards 

The  decimal  point  : 

ra.st  always  fail 

four  digits  from  the  end  of 

the  coefficient**  field 

The  decimal  toiut 

is  not  torched  in 

;  the  card;  the  computer  is 

programmed  to  provide  it  internally. 

Minus  signs  can  be 

punched  before  the  first  significant  digit 

but  >nder  no  circumstances 

can  minus  signs 

be  placed  outside  the  field 

occupied  by  the  coefficient.  Thus,  negative  coefficients  can  have  only 

three  significant  figures  before  the  decimal  point. 

It  cases  where  the  significant  digits  stop  before  the  end  of  the 
Held,  zeros  m.,st  be  fl  Lied  in  until  the  remainder  of  the  field  is  full 
of  digits.  No  zeros  are  required  before  the  first  significant  digit  in 
order  to  fill  the  field. 

Examples: 

a)  For  the  coefficient  -436.2’,  punch 

f-  .3  .6. a. 3 .0 .0  I 

b)  For  the  coefficient  0.0025,  punch 

I - 1  4  I - 1—| _ _ i  p  I 
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c)  For  the  coefficient  3261.,  punch 


1  3  1  2  ,  6  ,  1  ! _ 1 _ 1 _ 1 — | 

d)  Fcr  the  coefficient  0.,  punch 

I _ 1 _ 1 - 1 - 1 - i  i  J - 1 

3.4.6.  Order  of  the  Cards 

The  data  cards  are  placed,  in  the  order  described  above,  behind 
the  Fortran  program  deck  and  processed  as  a  standard  Fortran  compilation/ 
computation  sequence  with  either  the  1401  Fortran  compilation  deck  or 
1401  Fortran  tape. 

з. 5.  The  Output  Format 

Append!  x  II  shows  a  sequence  of  solution  sheets  to  a  typical  example 
of  e  quadratic  program  used  to  maximize  profits  ir.  a  firm.  The  first  sheet 
show*  the  initifti  data  ir.  the  solution  matrix  ar.l  the  initial  X  vector. 

The  l**t  sheet  shows  the  feasible  optima*  solution  zx.d  the  final  solution 
matrix 

By  suitable  m .difi-atior.*  to  toe  Joxtxaa  program,  the  feasible 
optimal  solution  may  he  placed  or  tap**  to  t«  -sed  it  furtler  processing 
by  another  program  («uzl  as  a  booickzepitg  routine)*  Baas,  the  quadratic 
program  may  fee  used  as  a  "subroutine*  fox  a  more  iaeolmed  probuer.. 

и.  Summary 

a)  'Tne  optimal  use  of  limited  resources,  where  constant  returns  to 
scale  cannot  be  assumed,  is  a  major  concern  in  the  theory  of  the  fires. 

b)  Quadratic  programming  is  a  practical  answer  to  this  problem. 
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c)  Medium-sized  digital  computers  are  available  to  a  great  many  firms. 
By  a  modification  of  the  Dantzig  quadratic  algorithm,  these  machines  may  be 
used  to  solve  a  great  many  of  the  optimization  problems  in  business  firms. 

d)  This  paper  describes  the  modified  Dantzig  quadratic  program. 

e)  A  Fortran  digital  computer  program  is  included,  with  instructions 
for  its  use,  and  a  practical  example. 
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Appendix  I 


PARAMI9R0206PS 

READ  5,  NT,  MT 
JB=NT+1 
M=NT+MT 
N=M+M 

DIMENSION  A(30,15),C(15),X(30) 
PRINT  8 
DO  3  1=1,30 
X(l)=0. 

DO  2  J=l,15 

2  A(l,J)=0. 

3  CONTINUE 

READ  6,(C(J),J=1,10) 

READ  7,(C(J),J=11,15) 

DO  5  1=1, NT 

5  READ  6,(A(I,J),J=1,10) 

DO  60  J=JB,M 
READ  6, (A(l ,J), 1=1, 10) 

X(J)=C(J) 

MR=J+M 
DO  50  1=1, NT 
50  A(MR,I)=A(I,J) 

60  CONTINUE 
DO  70  J=1,NT 
l=M+J 
X(I)=-C(J) 

70  A(I,J)=-1. 

PRINT  600, (CO), 1=1 ,15) 

DO  200  1=1, N 

200  PRINT  615 , (A(l ,J),J=1 ,15),X(I) 

C  FIND  INCOMING  ACTIVITY,  K 
JB=0 


1-1 


MR=M+1 
105  E=0. 

K=0 

DO  120  l=MR,N 
IF(E-Xd))  120,120,110 
110  E=X(I) 

K=l-M 

120  CONTINUE 
JA=M+K 

C  OPTIMALITY  CHECK 

IF(X(JA))  205,700,700 
C  FIND  BINDING  CONSTRAINT,  L 
205  L=0 

£=999999. 

DO  230  J=1,M 

IF(A(JA,J)*C(J))  219,204,219 
204  IF(X(J))  219,230,219 

219  IF(A(K,J))  212,230,209 
212  IF(C(J))  210,230,230 

209  IF(C(J))  230,230,210 

210  D=C(J)/A(K,J) 

IF(E-D)  230,230,220 

220  E=D 
L=J 

230  CONTINUE 

IF(L)  235,235,240 
C  FIND  DROPPING  ACTIVITY 

240  DO  243  1=1, N 
IF(X(D)  241,243,241 

241  IF(A(I,L))  242,243,242 

242  KA=I 

243  CONTINUE 
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MODIFY  A  MATRIX  AND  SOLVE  FOR  NEW  CONSTRAINTS,  C(J) 
E=A(K,L) 

DO  300  l=l,N 
300  A(I,L)=A(I,L)/E 
C(L)=C(L)/E 
DO  320  J=1,M 
IF(J-L)  310,320,310 
310  D=A(K,J> 

C(J)=C(J)-D*C(L> 

DO  315  1=1, N 
315  A(l,J)=A(l,J)-D*A(l,L> 

320  CONTINUE 
JB=JB+1 

PRINT  630,JB,K,KA,L 
PRINT  600,(C(l),l=l,15) 

C  SOLVE  FOR  NEW  ACTIVITY  LEVELS 

X(K)=C(L) 

X(KA)=0. 

DO  400  1=1, N 
IF(X(D)  330,400,330 
330  E=0. 

DO  340  J=1,M 
340  E=E+C(J)*A(I,J) 

X(I)=E 

400  PRINT  615  ,(A(I  ,J),J=1 ,15),X(I) 

C  COMPLEMENTARITY  CHECK 

IF(X(JA»  510,595,595 
510  K=KA+M 
GO  TO  205 
595  PRINT  640 
GO  TO  105 
700  PRINT  710 
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STOP  222 
235  STOP  888 

5  FORMAT  (212) 

6  FORMAT  (10F8.4) 

7  FORMAT  (5F8.4) 

8  FORMAT  (32H1WOLFE  SIMPLEX  QUADRATIC  PROGRAM,//) 

600  FORMAT  <1X,15F8.4,//) 

615  FORMAT  (1X,16F8.4) 

630  FORMAT  (12H1IN  SOLUTION ,14,  9H,ACTIVITY,I4,6HIS  IN,,I4,21HDR0PS,  AN 
ID  CONSTRAINT, 14 ,6HBINDS.) 

640  FORMAT  (28H0THIS  IS  A  FEASIBLE  SOLUTION) 

710  FORMAT  (21H  AND  OPTIMAL  SOLUTION) 

END 
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